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Abstract 

Within the real-time formulation of non-equilibrium field theory generalized trans- 
port equations are derived avoiding the standard quasiparticle approximation. They 
permit to include unstable particles into the transport scheme. In order to achieve 
a self-consistent, conserving and thermodynamically consistent description, we gen- 
eralize the Baym's ^-functional method to genuine non-equilibrium processes. This 
scheme may be closed at any desired loop order of the diagrams of the functional 
<3? this way defining a consistent effective theory. By means of a first-order gradient 
approximation the corresponding Kadanoff-Baym equations are converted into a 
set of coupled equations. This set consists of a time-irreversible generalized kinetic 
equation for the slowly varying space-time part of the phase-space distributions and 
a retarded equation, which provides the fast micro-scale dynamics represented by 
the four-momentum part of the distributions. Thereby, no constraint to the mass 
shell of the particles is required any further and the corresponding spectral mass 
distributions are treated dynamically. The description naturally includes all those 
quantum features already inherent in the corresponding equilibrium limit (Mat- 
subara formalism). Memory effects appearing in collision term diagrams of higher 
order are discussed. The variational properties of ^-functional permit to derive a 
generalized expression for the non-equilibrium kinetic entropy flow, which includes 
corrections from fluctuations and mass width effects. In special cases an //-theorem 
can be demonstrated implying that the entropy can only increase with time. Mem- 
ory effects in the kinetic terms provide corrections to the kinetic entropy flow that in 
equilibrium limit recover the famous bosonic type T 3 In T correction to the specific 
heat of Fermi liquids like Helium-3. 



1 Introduction 



A proper dynamical scheme in terms of a transport concept that deals with unstable 
particles, like resonances, i.e. particles with a broad mass width, is still lacking. Rather ad- 
hoc recipes are in use that sometimes violate basic requirements as given by fundamental 
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symmetries and conservation laws, detailed balance or thermodynamic consistency. The 
problem of conserving approximation has first been addressed by Baym and Kadanoff 
[1,2]. They started from an equilibrium description in the imaginary time formalism and 
discussed the response to external disturbances. Baym, in particular, showed [2] that 
any approximation, in order to be conserving, must be so-called ^-derivable. It turned 
out that the <£> functional required is precisely the auxiliary functional introduced by 
Luttinger and Ward [3] (see also ref. [4]) in connection with the thermodynamic potential. 
In nonequilibrium case the problem of conserving approximations is even more severe than 
in near-equilibrium linear-response theory. 

The appropriate frame for the description of nonequilibrium processes is the real-time 
formalism of quantum field theory, developed by Schwinger, Kadanoff, Baym and Keldysh 
[5-7]. In ref. [8] we have introduced the generating functional $ on the real-time contour. 
It is determined by the sum of all closed vacuum skeleton diagrams in terms of full classical 
fields and full Green's functions. All important quantities of a system (such as the sources 
of classical fields, self-energies, interaction energy, etc.) are derived by variations of the $ 
functional with respect to its arguments. The advantage of the $ functional is that one 
may formulate various approximations in terms of approximate $ (so called ^-derivable 
approximations), which preserve the conservation laws related to global symmetries of the 
underlying theory and thermodynamic consistency. Thereby, one may restrict oneself to 
either few diagrams only or to some sub-set of diagrams for $. Note that far not every 
approximation scheme, e.g., at the level of self-energies, possesses such properties. 

Functional variational methods allow to derive the corresponding Dyson equations on real- 
time contour. The steps towards generalized kinetic equations are provided by the Wigner 
transformation of all two-point functions and by the first-order gradient approximation. 
The former formulates all quantities in terms of phase-space distribution functions in four 
dimensions, i.e. as a function of energy and momentum for any space-time coordinate. 
The gradient approximation leads to Poisson-bracket expressions which permit a classical 
interpretation. For complicated collision terms, the latter also give rise to memory effects. 
The quasiparticle approximation is not required at any step. Rather one can keep the finite 
damping width of the particles. The corresponding dynamical information contained in the 
spectral functions is determined by local retarded equations. We derive closed expressions 
for the nonequilibrium energy momentum tensor in the ^-derivable scheme, which also 
include terms arising from the fluctuations in the system and the finite mass width of the 
particles. 

The second main issue of this paper is the role of entropy. Although the entropy is a central 
quantity in thermodynamics and statistical mechanics, many problems connected with it, 
in particular, its description in terms of Green's functions in the nonequilibrium case, are 
still open. One can find related discussions in many textbooks and reviews, e.g., in refs 
[9-13]. The thermodynamic entropy has extensively been discussed in the literature at the 
end of sixties and beginning of seventies, cf. refs [14-16] and references therein. The fact 
that the generalized kinetic equation possesses a proper thermodynamic limit, does not 
yet imply that this limit will be approached during the evolution. The latter is however 
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ensured, if one can prove an if -theorem for the equations of motion. Using the ^-derivable 
properties we are able to get an expression, which takes the sense of a nonequilibrium 
kinetic entropy expressed in terms of Green's functions and self-energies. It generalizes to 
nonequilibrium the well known equilibrium expression obtained in Matsubara formalism 
within the ^-functional scheme. We derived conditions for an if -theorem, which apply in 
many cases discussed in transport problems. In addition, the kinetic entropy expression 
is extended beyond the quasiparticle limit thus accounting for the fluctuations. Memory 
effects lead to further corrections to the entropy flow. The latter have to be established 
for each particular case. In spite of many attempts so far, if -theorem has not yet been 
proven even for the classical kinetic equation including triple collision term, cf. ref. [13]. 
Therefore, up to now there exists no appropriate kinetic entropy expression (derived from 
generalized kinetic description). Our demonstration of if -theorem for particular cases 
including the fluctuation parts may be viewed as step towards this general goal. 

The paper is organized as follows. Within the real-time contour formulation of nonequi- 
librium dynamics we define a generating functional T together with the auxiliary $ func- 
tional of Baym, cf. ref. [8], (sect. 2). Then, for slowly varying generalized distribution 
functions we derive a generalized kinetic equation consistent with $, using the gradient 
approximation (sect. 3). Arguments are given that the generalized distribution functions 
remain positive during the time evolution and thus permit a probabilistic interpretation 
in four-momentum space at each space-time point. The properties of <£> determine the 
conservation laws in terms of conserved Noether currents and energy-momentum tensor. 
Further properties are exploited with the help of the decomposition rules for the diagrams 
of $ formulated in terms of full Green's functions [17]. These rules lead to a multi- process 
decomposition of a ^-derivable collision term, which contains space-time local and non- 
local (memory) parts (sect. 4). Time-irreversibility is discussed in connection with our 
generalized kinetic description, and expression for the kinetic entropy flow is derived. In 
specific cases an explicit growth of the kinetic entropy with time, i.e. an if -theorem (sect. 
5) is demonstrated. The proof once again relies on the $-derivability of the collision term. 
Markovian and memory contributions to the entropy flow are discussed at a specific ex- 
ample. Finally, the proper thermodynamic limit of the kinetic entropy is demonstrated 
(sect. 6). Some formal details are deferred to Appendices. 

To be specific, we concentrate on systems of nonrelativistic particles. Bosonic mean fields 
are not treated in this paper; they can however be included along the lines given in ref. 
[8]. Whenever possible we omit the index a denoting the different particle species and 
intrinsic quantum numbers in a multi-component system. 
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2 Generating Functionals T and $ 



2.1 $ Functional and Dyson's Equation on Real- Time Contour 

We assume the nonequilibrium system to be prepared at some initial time to i n terms 
of a given density operator p — J2 a Pa \&) (&\, where the |a) form a complete set of 
eigenstates of po- 

All observables can be expressed through n-point Wightman functions of Heisenberg op- 
erators A(ti), . . . , 0{t n ) at some later times 

(6(t n ) . . . B(t 2 )A(ti)^ =: Tr 6(t n ) . . . B^A^p^h) 

= Y,Pa (a\ 0(t n ) . . .Bit^Ait,) \a) . (2.1) 

a 

Note the fixed operator ordering for Wightman functions. 




to it oo 



Figure 1: Closed real-time contour with two external points x,y on the contour. 

The nonequilibrium theory can entirely be formulated on a single special contour, the 
closed real-time contour (see figure 1) with the time argument running from t to oo along 
time-ordered branch and back to to along anti-time- ordered branch. The contour Green's 
function is defined as the expectation value of contour-ordered products of operators 

iG(x,y) = {T c <p(x)<pi(y)), (2.2) 

where Tq denotes the special time-ordering operator, which orders the operators according 
to a time parameter running along the time contour C. With the aim to come to a self- 
consistent treatment in terms of two-point Green's functions, one has to ignore higher 
order correlations which would be present in the exact dynamical equation for two-point 
Green's functions. This step is achieve by the Wick decomposition which leads to the 
Dyson equation now formulated on the real time contour 

S x G(x, y) = 5 c (x, y) + J dzS(x, z)G(z, y), (2.3) 

c 
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S x G°(x,y) = 5 c (x,y), 
where 



(2.4) 



S. = ift + ^fl£ (2.5) 

in nonrelativistic kinematics. Here 5c{x,y) is ^-function on the contour, G° is the free 
Green's function and S is the self-energy, expressed via one-particle Green's function in 
the usual way. This way one ignores higher order correlation functions, cf. refs [18,19]. 
Then the self-energy itself becomes a functional of the Green's functions (2.2) with the 
consequence that the standard Wick decomposition can be applied for the self-energy 
and one thus obtains standard contour diagrams, cf. ref. [8]. In such a diagrammatic 
representation eq. (2.3) is given by 




+ — - — f-iEJ — «— (2.6) 



with the two-point functions — iS {G(x,y)}, as driving terms, and the Green's functions 



iG°(x,y)= — - — , iG(x,y)= — . (2.7) 

v ' y; x y ' v yj x y v ' 

The arrow always points towards the <^-field operator in the contour ordered expression, 
cf. eq. (2.2). 

The physical condition, which justifies to disregard higher order correlations, is that the 
typical interaction time T mt for the change of the correlation functions is much less than 
the typical relaxation time r re i, which determines the system evolution. Describing the 
system at times t — to ^> T int , one can neglect initial correlations, since they are supposed 
to be dying out at times ~ r int in accordance with the Bogolyubov's principle of the 
weakening of initial correlations.^ This coarse-graining leads to time-irreversibility, as we 
will see below. Alternatively, one also could suppose the initial state to be uncorrelated, 
like an equilibrium ideal gas state, which corresponds to an information loss right from 
the beginning, cf. refs [6,20]. 

In the imaginary-time formalism, Luttinger and Ward [3] have shown that the generating 
functional T can be expressed in terms of an auxiliary functional $, where $ is solely 
given in terms of full, i.e. resummed, propagators. Following ref. [8], it can be generalized 
to the real-time case with the following diagrammatic representation 



1 Actually, considering dilute systems, Bogolyubov suggested the weakening of all correlations, 
whereas we use a weaker assumption on the weakening of only short-time (~ r; n t) correlations, 
cf. ref. [12]. 
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ir{G,A} = ir°{G' } 




± In (l - OG E) ±0G0S + i${G,A} 

Here upper signs relate to fermion quantities, whereas lower signs, to boson ones. For the 
sake of mathematical convenience, cf. ref. [8], a scaling parameter A is introduced in each 
interaction vertex. For the physical case, A = 1. The integer ns counts the number of self- 
energy insertions in the ring diagrams. For the closed diagrams of <3> the value ri\ counts 
the number of vertices building up each such diagram. Due to the global factor l/ri\ the 
set of ^-diagrams is not resumable in the standard diagrammatic sense. For nonrelativistic 
instantaneous two-body interaction each interaction counts as two vertices, cf. Appendix 
C. The term T solely depends on the free propagator G°. The diagrams contributing to 
$ are given in terms of full propagators G (thick lines). These diagrams are two-particle 
irreducible (label c2), i.e. they cannot be decomposed into two pieces by cutting two 
propagator lines. The latter property matches diagrammatic rules for the resummed self- 
energy E(x, y), which results from functional variation of <3> with respect to the propagator 
G(y,x), i.e. 



- i£ = zf5i<S>/5iG. (2.9) 

Relation (2.9) directly follows from the stationarity condition of T (2.8) with respect to 
variations in the full propagator G on the contour 



5T {G, A} /SG = 0, (2.10) 

which also provides the Dyson's equation (2.3). In graphical terms, the variation (2.9) 
is realized by opening a propagator line in any diagram of $. The resulting set of thus 
opened diagrams must then be that of proper skeleton diagrams of £ in terms of full 
propagators, i.e. void of any self-energy insertions. 

Subsequently we will not use the contour representation but rather its decomposition into 
the two branches the time-ordered and the anti-time ordered branch using a { — h} matrix 
notation as explained in Appendix A, cf. also refs [21,8]. 
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2.2 $ -Derivable Approximation Scheme 



For any practical calculation one has to truncate the scheme. In the weak-coupling limit, 
the perturbative expansion may be restricted to a certain order. Then no particular prob- 
lems are encountered as far as conservation laws are concerned, since they are fulfilled 
order by order in perturbation theory. On the other hand, such perturbative expansion 
may not be adequate, as, for example, in the strong coupling limit, where resummation 
concepts have to be applied. Such schemes sum up certain sub-series of diagrams to any 
order. Furthermore, with the aim to solve dynamical equations of motion, such as trans- 
port equations, one automatically resums all terms in the equations of motion to any 
order. For such resumming schemes, the situation with conservation laws is not quite as 
clear. 

We consider the so-called ^-derivable approximations, first introduced by Baym [2] within 
the imaginary time formulation. A ^-derivable approximation is constructed by confining 
the infinite set of diagrams for <3> to either only a few of them or some sub-series of them. 
Note that $ itself is constructed in terms of "full" Green's functions, where "full" now 
takes the sense of solving self-consistently the Dyson's equation with the driving terms 
derived from this $ through relation (2.9). It means that even restricting ourselves to a 
single diagram in <3>, in fact, we deal with a whole sub- series of diagrams in terms of free 
Green's functions, and "full" takes the sense of the sum of this whole sub-series. Thus, 
a ^-derivable approximation offers a natural way of introducing closed, i.e. consistent 
approximation schemes based on summation of diagrammatic sub-series. In order to pre- 
serve the original symmetry of the exact $, we postulate that the set of diagrams defining 
the ^-derivable approximation complies with all such symmetries. As a consequence, ap- 
proximate forms of $( appr -) define effective theories, where <|>( appr -) serves as a generating 
functional for approximate self-energies £( appr -)(a;, y) 



which then enter as driving terms for the Dyson's equations (2.3). The propagators, solving 
this set of Dyson's equations, are still called "full" in the sense of the $( appr -)-derivable 
scheme. Below, we omit the superscript "appr." . 



3 Generalized Kinetic Equation 

3.1 Gradient Expansion Scheme 

For slightly inhomogeneous and slowly evolving systems, the degrees of freedom can be 
subdivided into rapid and slow ones. Any kinetic approximation is essentially based on 



£i$(appr.) 



(2.11) 
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this assumption. Then for any two-point function F(x,y), one separates the variable 
£ = (£1 — t 2 , t*i — 7*2), which relates to rapid and short-ranged microscopic processes, and 
the variable X = |(ti+t 2 , r 1 +r 2 ), which refers to slow and long-ranged collective motions. 
The Wigner transformation, i.e. the Fourier transformation in four-space difference £ = 
x — y to four-momentum p leads to the corresponding Wigner densities in four-phase- 
space. Since the Wigner transformation is defined for physical space-time coordinates 
rather than for contour coordinates one has to decompose the contour integrations into 
its two branches, the time-ordered { — } branch and the anti-time ordered {+} branch. 
This is explained in detail in Appendix A. Two-point functions then become matrices of 
the contour decomposed { — h} components with physical space-time arguments. Thus 

FV(X;p) = J d£ e i ^(X + £/2,X-£/2), i,j G {-+} (3.1) 



leads to a four-phase-space representation of two-point functions. The retarded or ad- 
vanced relations between two-point functions formulated on real-time contour and those 
in matrix notation, F 1 *, presented in Appendix A, are thereby preserved. 

The Wigner transformation of Dyson's equation (2.3) is straight forward. Taking the 
difference and half-sum of Dyson's equation (2.3) and the corresponding adjoint equation 
and doing the Wigner transformation we arrive at equations 



iv li d$G i *(X,p) = J d&rtj dz(E(x i ,z)G(z,y^-G(x i ,z)E(z,y^) , (3.2) 

c 

Q x GV(X,p) = aV + ±fdZt**f dz(E(x i ,z)G(z,y^+G(x i ,z)E(z,y^) , (3.3) 

c 

where <r u accounts for the integration sense on the two contour branches, cf. eqs (A. 2), 
(A. 3). For nonrelativistic kinematics = (l,p/m), and Qx = Po — p 2 /2m — d\/d>m. 
In this matrix notation, two of equations (3.2) and (3.3), involving G~ + and G + ~ on 
the left-hand side, are known as Kadanoff-Baym equations in Wigner representation [6]. 
Particular combinations of these equations lead to the retarded and advanced equations 
which completely decouple and involve only integrations over physical times rather than 
contour times. 

Standard transport descriptions usually involve two further approximation steps: (i) the 
gradient expansion for the slow degrees of freedom, as well as (ii) the quasiparticle approx- 
imation for rapid ones. We intend to avoid the latter approximation and will solely deal 
with the gradient approximation for slow collective motions by performing the gradient 
expansion of eqs (3.2) and (3.3). This step preserves all the invariances of the <E> functional 
in a ^-derivable approximation. 
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Technically, the gradient expansion is performed as follows. The Wigner transformation 
of a convolution of two two-point functions is given by 



d£e 



dzf(x,z)<p(z,y) 



cxp 



ill 



(d p d x > - d x d p/ ) 



f(X, P )<p(X', P ') 



pi=p,X'=X 



Here, we have temporally restored h for the sake of physical lucidity Expanding this 
expression in powers of hdx • d p to the first order, we obtain 

J d£e«(/ dzf(x,zMz,y)^ ~ f(X,pMX,p) + ^{f(X,p),< P (X,p)}, (3.5) 



where 

{f(X,p)MX,p)} = - (3.6) 

is the classical Poisson bracket in covariant notation. We would like to stress that the 
smallness of the hdx • d p comes solely from the smallness of space-time gradients dx, 
while momentum derivatives d p are not assumed to be small! They must be not too large 
only in order to preserve the smallness of the total parameter hdx • d p . This point is 
sometimes incorrectly treated in the literature. 



3.2 Kinetic Equation in Physical Notation 



It is helpful to avoid all the imaginary factors inherent in the standard Green's function 
formulation and introduce two quantities which are real and, in the quasi-homogeneous 
limit, positive, cf. subsect. 3.5, and therefore have a straightforward physical interpreta- 
tion, much like for the Boltzmann equation. We define 

F(X,p) = A(X,p)f(X,p) = ( T )iG- + (X,p), 

F(X,p) = A(X,p)[l T f(X,p)\ = iG + -(X,p), (3.7) 

for the generalized Wigner functions F and F and the corresponding /onr-phase-space 
distribution functions f(X,p) and Fermi/Bose factors [1 =F f(X,p)}. Here 



A(X, p) = -21mG R {X, p)=F±F = i (g + ~ - CT+) (3.8) 

is the spectral function, where G R is the retarded propagator, cf. eq. (A. 6) of Appendix 
A. According to relations (A. 6) and (A. 7) between Green's functions G lj , only two real 
functions of all these G 1 ^ are required for a complete description of the system 's evolution. 
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The reduced gain and loss rates of the collision integral are defined as 



T in (X,p) = T(X,p) 1 (X,p) = ( T )iE- + pT,p), (3.9) 
T out (X,p) = T(X,p)[lTl(X,p)] =i£ + -pT,p) (3.10) 

with the damping width 

T(x,p) = -2imi: R (x,p) = r out (x, P )±r in (x, P ), (3.11) 

where S K is the retarded self-energy, cf. eq. (A. 6). The opposite combination 

i(x, P ) = [2 7 t i]r = r in (x, P ) =f r out (x, P ), (3.12) 

is related to fluctuations. The dimensionless quantity 7 is introduced for further conve- 
nience. 

In terms of the new notation (3.7)— (3.11) and within the first-order gradient approxima- 
tion, the equations (3.2) for F and F take the form 

VF(X,p) - {T in ,ReG R } = C(X,p), (3.13) 
VF(X,p) - {T out ,ReG R } = TC(X,p) (3.14) 

which we denote as the transport equations. Here the differential drift operator is defined 

as 



Furthermore 



cpr, P ) = r in (x,p)F(x, P ) - r out (x,p)F(x,p) = Ar[ 7 - /] (3.16) 

is the collision term with the dimensionless functions / and 7 defined in (3.7) and (3.9), 
while jr in , ReCr^j and |r out ,ReG fi | are the fluctuation terms. Within the same approx- 
imation level eq. (3.3) provides us with two alternative equations for F and F 



MF - ReG R r in = l - ({r, F} - {T in , A}) , (3.17) 
MF - ReG R T out = - ({r, P) - {r out , A}) (3.18) 
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with the "mass" function 



M(X,p) = p - ^P 2 - ReE R (X,p), (3.19) 

which relates to the drift operator via Df = {M, /} for any four-phase-space function 
/. Eqs (3.17), (3.18) are reported as the mass shell equations. Appropriate combinations 
of the two sets of equations (3.13)-(3.14) and (3.17)-(3.18) provide us with the retarded 
equations 



VG R (X,p)+ 1 - {T,G R } = 0, (3.20) 
(M(X,p) + ^(X,p)) G R (X,p) = 1. (3.21) 
This subset (3.20)-(3.21) is solved by [22] 



G R = 



M(X,p)+iT(X : p)/2 



A(X,p) 



ReG R (X,p) 



M 2 (X,p) + r 2 (X,p)/4' 

M(X,p) 
M 2 (X,p) + T 2 (X,p)/4 



(3.22) 



The spectral function satisfies the sum-rule 



J^A(X,p) = l, (3.23) 

— oo 

which follows from the canonical equal-time (anti) commutation relations for (fermionic) 
bosonic field operators. 

With the solution (3.22) for G R equations (3.13) and (3.17) b ecome identical to (3.14) 
and (3.18), respectively. However, equations (3.13) and (3.17) are nonidentical, while they 
were identical before the gradient expansion. We will now discuss under which conditions 
they are equivalent. First, this takes place in the standard quasiparticle limit when one 
puts r — > in the Green's functions. This limit will be discussed in the next subsect. 
The equivalence is also given in the case we are interested in here, where T is of finite 
value, however, one needs \f — j\ <C 1 which does not necessarily require the system to 
be close to local equilibrium. In both cases the gradient approximation implies that the 
macroscopic time scale r macro , characterizing kinetic processes, is much larger than the 
microscopic time scale r m j cro , relating to rapid microscopic processes. 

Using the limit \f(X,p) — ^y(X,p)\ <^ 1 and writing all r in = 'yT terms as (/ + (7 — f))T, 
the transport equation (3.13) and the mass-shell equation (3.17) take the forms 
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y (T {M, f} — M {r, /}) = TA( 7 - /) - r kin , (3.24) 

— ^— (r {M, /} - M {r, /}) = AM( 7 - /) - r mass _ eq , (3.25) 

where the remaining terms r kin and r mass _ shell are of order (7 — /) times gradient terms, 
i.e. of second order in the gradient expansionP]. From eq. (3.22) one sees that TR,eG R = 
MA. Using this we observe that the mass-shell equation loses its original sense, since to 
leading order relation (3.25) is equivalent to the kinetic equation (3.24). However, the still 
remaining difference in the second-order terms is inconvenient from the practical point of 
view. Following Botermans and Malfliet [22] who first suggested to drop the above r^ in 
term in the kinetic equation we put r^n = r mass = 0, since then both equations indeed 
become identical. This step amounts to replace the r in and r out terms by fT and (1 =F/)r 
in all Poisson brackets. The so obtained generalized kinetic equations for F and F then 
read 



VF(X,p) - |r^,ReG R } = C(X,p), (3.27) 
VF(X,p) - S^j,ReG R \=TC(X,p). (3.28) 

Although these two equations are identical to each other, we nevertheless have presented 
both of them for the sake of further convenience. With the help of the retarded relations 
(3.22) both equations can be converted to one generalized transport equation for the 
phase-space occupation functions / 



^-(TVf-M{TJ}) = C. (3.29) 
To get eq. (3.29) we used that the additional Poisson-bracket fluctuation term becomes 

{r/.Rec-} = v(T f) + (M f; r y /4)2 {r, /} . 



We now provide a physical interpretation of various terms in the generalized kinetic equa- 
tion (3.27) or (3.29). This physical interpretation relies on the similarity of most of the 
terms to conventional kinetic equations, for example, such as the Landau equation for 

2 These terms have the explicit form 

r k in = - {F(7 - f),ReG R } , r mass _ eq = J {r( 7 -f),A}. (3.26) 
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Fermi liquids (see, e.g., refs [16,21]) and ref. [23] for the relativistic version), proposed by 
Landau on the basis of an intuitive quasiparticle picture [24]. 

For this purpose it is advantage to convert the drift operator (3.15) into a space and time 
separated form 



takes the meaning of the group 3-velocity in the quasiparticle approximation, while Z is 
the standard renormalization factor. 

Thus, the drift term T>F on the l.h.s. of eq. (3.27) is the usual kinetic drift term including 
the corrections from the self-consistent field Re£ R into the convective transfer of real and 
also virtual particles. In the collisionless case VF = (Vlasov equation), the quasi-linear 
first-order differential operator V defines characteristic curves. They are the standard 
classical paths. The four-phase-space probability F(X,p) is conserved along these paths. 
The formulation in terms of a Poisson bracket in four dimensions implies a generalized Li- 
ouville theorem. In the collisional case, both the collision, C, and fluctuation terms (3.30) 
change the phase-space probabilities of the "generalized" particles during their propaga- 
tion along the "generalized" classical paths given by V. We use the term "generalized" in 
order to emphasize that particles are no longer bound to their mass-shell, M — 0, during 
propagation due to the collision term, i.e. due to decay, creation or scattering processes. 

The r.h.s. of eq. (3.27) specifies the collision term C in terms of gain and loss terms, 
which also can account for multi-particle processes. Since F includes a factor A, C fur- 
ther deviates from the standard Boltzmann-type form in as much that it is multiplied by 
the spectral function A, which accounts for the finite width of the particles. The Poisson- 
bracket term (3.30) is special. It contains genuine contributions from the finite mass width 
of the particles and describes the response of the surrounding matter due to fluctuations. 
This can be seen from the conservation laws discussed below. In particular the first term 
in (3.30) gives rise to a back-flow component of the surrounding matter. It restores the 
Noether currents as the conserved ones rather than the intuitively expected sum of con- 
vective currents arising from the convective T>F terms in (3.27). In eq. (3.29) one further 
realizes that the function \A 2 Y in front of the drift term is more sharply peaked than the 
original spectral function A. Both \A 2 T and A are reduced to the same (^-function in the 
quasiparticle limit. The M {T, /} term, which gives no contribution in the quasiparticle 
limit due to the factor M, represents a specific off-mass-shell response, cf. ref. [25]. 

So far the gradient approximation has been applied to the space-time foldings occurring 



V = % (dt + v g dx) + d t ReZ R ■ d Po - dxReE^ • dp , 



(3.31) 



where 




with 




(3.32) 
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between the self-energies and the propagators appearing in the collision term. They give 
rise to the gradient terms on the l.h.s. of the kinetic equations. This is sufficient as long 
as the self-energies are calculated without further approximation. Commonly one likes 
to obtain also T in and r out in a kind of local approximation evaluated with all Green's 
functions taken at the same space-time point X. Then for self-energy diagrams with more 
than two points, also nonlocal gradient corrections within these diagrams arise, which 
have to be accounted for in a consistent gradient approximation scheme. The latter fact 
gives rise to memory effects, which will be discussed in detail in subsect. 4.2, below. 

In all derivations above the r in and r out terms in all Poisson brackets have consistently 
been replaced by fT and (l^ff)T, respectively. In fact, this procedure is optional. However, 
in doing so we gained few advantages. First, all our generalized kinetic equations, (3.24), 
and mass-shell equations, (3.25) are exactly equivalent to each other, as they were in exact 
theory before the gradient expansion. Second, the so obtained generalized kinetic equation 
(3.29) has particular feature with respect to the definition of a nonequilibrium entropy 
flow in connection with the formulation of an exact H-theorem in certain cases (sect. 5). 
If we refrained from these substitutions, these both features would become approximate 
with the accuracy up to second-order gradients. While this accuracy still complies with 
the gradient approximation level, it makes our first-order approximation less elegant. 

3.3 Quasiparticle Limit 

The quasiparticle picture is recovered in the limit of small (while nonzero) width V . Small 
means that Y should be less than all typical scales in the problem. This allows to put 
T zero in all Green's functions which are, thereby, called quasiparticle Green's functions, 
nonzero means that V is kept finite in the collision term being now calculated with the 
help of these quasiparticle Green's functions. 

Thus, in this limit the spectral function (3.22) takes the form 

A qp = 2ir signr 5{M) (3.33) 

with M given by eq. (3.19). As the right hand side of the mass shell equation (3.17) 
contains one V factor more than the left hand side, it reduces to MA(f — 7) = requiring 
M to vanish within the range of the spectral function. This just complies with the on- 
shell 5-function in eq. (3.33) which with M = provides the dispersion equation for the 
quasiparticle energies e(X,p) 

e(X,p) = ^-p 2 + Re S fi (X,e(X,p),p) . (3.34) 

This dispersion equation may have several solutions (branches). For brevity we drop a 
possible branch index. In the quasiparticle limit the generalized Wigner functions and the 
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Fermi/Bose factors (3.7) are presented in the form 



F^(X,p) = 2vr Z*>(X,p) 5(p - e(X,p)) / qp (X,p), 
F»(X,p) = 2tt Z™(X,p) 5 (p - e(X,p)) [1 T f qp (X,p)] , 



(3.35) 
(3.36) 



with the quasiparticle distribution functions 



r*(X,p) = f(X,e(X,p),p), 



(3.37) 



which depend on the 3-momentum only rather than on the 4-momentum, while 



denotes the quasiparticle normalization factor (cf. eq. (3.32)). In practice, Z qp is always 
positive for truly nonrelativistic systems [26]. We have replaced p by e(X,p) in all rela- 
tions above due to the on-shell 5-function in eq. (3.33). 

Integration of eq. (3.27) over p Q in a narrow region of width T around the branch position 
e(X,p), where one can identify A = A qp , leads to the standard quasiparticle transport 
equation 



dtfP(X,p) + [ r P) d^(X,p) - dxe(X,p) df l A ' P) = C qp (X,p). (3.39) 



Here C qp (X, p) is the corresponding collision integral now expressed in terms of quasi- 
particle distributions. The fluctuation term jr/, ReG^j integrated over p becomes of 
second order in T and thus drops out of this quasiparticle equation. 

This quasiparticle transport equation is consistent as long as its evolution does not lead 
beyond its region of validity, i.e. into region where T becomes comparable to typical energy 
scales in the problem. For instance, the consistency of the Landau Fermi-liquid theory 
is preserved, provided only low-excited Fermi systems are considered. Then the Pauli- 
blocking prevents particles to scatter to states laying far away from the Fermi surface, i.e. 
to a region where the quasiparticle approximation would be violated. 

In above approach neither any uncontrolled approximations nor any specific ansatz has 
been introduced. The quasiparticle limit of the generalized kinetic equation directly follows 
from our consideration and encounters no problems for the eliminating the extra Poisson 
bracket term. Here a caution should be expressed. The quasiparticle form of the spectral 
function of eq. (3.33) gives 



Z qp (X,p) = Z(X,e(X,p),p) 



(3.38) 






(3.40) 



— oo 
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for the sum- rule rather than unity. The Z qp factor appears here, because in (3.40) we 
artificially extended the quasiparticle form of A to all energies and did not account for 
background terms outside the quasiparticle region. However, in this form the sum rule 
is often taken as the canonical one for quasiparticles, cf. ref. [27], in order to construct 
an effective theory only in terms of quasiparticles. This recipe works, if the quasiparticle 
energy region is indeed well separated from the region, where particles have a rather large 
width. 



3.4 Conservations of Charge and Energy-Momentum 



Transport equation (3.27) weighted either with the charge e or with 4-momentum p", inte- 
grated over momentum and summed over internal degrees of freedom like spin (Tr) gives 
rise to the charge or energy-momentum conservation laws, respectively, with the Noether 
4-current and Noether energy-momentum tensor defined by the following expressions 



nX)=eTrf-^- 4 v»F(X,p), 



9^(X) = Tr J 0- 4 v^F(X,p)+g^ (S^(X) - S pot (X)) 



(3.41) 
(3.42) 



Here 



(3.43) 

A=l 

is interaction energy density, which in terms of $ is given by a functional variation with 
respect to a space-time dependent coupling strength of interaction part of the Lagrangian 
density £ mt — > X(x)C mt , cf. ref. [8]. The potential energy density £ pot takes the form 



r-Hx) = <-2-po) = — 



£ P ot = ^ J 



d 4 p 



ReL R F + ReG R ^-F 
A 



(3.44) 



Whereas the first term in squared brackets complies with quasiparticle expectations, 
namely mean potential times density, the second term displays the role of fluctuations 
in the potential energy density. 

For specific interactions with the same number a of field operators attached to any vertex 
of £ mt , one simply deduces 



£™\X,p) = -£»°\X,p). 
a 



(3.45) 
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Please also notice from (3.42) that the special combination 



OO PO + ^E e "P0 = Tr/ ^4(X,p)/(X,p) 



(3.46) 



depends on the specific form of the interaction only via the spectral function. As we will 
see below, this combination simply relates to the entropy density in local thermodynamic 
equilibrium. 

The conservation laws only hold, if all the self-energies are ^-derivable. In ref. [8], it has 
been shown that this implies the following consistency relations obtained after Wigner 
transformation and first-order gradient expansion, (a) for the conserved current 



lTr /(2^[{ Re ^ F }-{ Re ^>} + C 



(2tt)< 



= 0, 



(3.47) 



and (b) for the energy-momentum tensor 



(3.48) 



The contributions from the Markovian collision term C drop out in both cases, cf. eq. 
(4.8) below. The first term in each of the two relations refers to the change from the free 
velocity v to the group velocity v g , cf. eq. (3.32), in the medium. It can therefore be 
associated with a corresponding drag-flow contribution of the surrounding matter to the 
current or energy-momentum flow. The second (fluctuation) term compensates the former 
contribution and can therefore be associated with a back-flow contribution, which restores 
the Noether expressions (3.41) and (3.42) to be indeed the conserved quantities. In this 
compensation we see the essential role of the fluctuation term (3.30) in the generalized 
kinetic equation. Dropping or approximating this term would spoil the conservation laws. 
Indeed, both expressions (3.41) and (3.42) comply exactly with the generalized kinetic 
equation (3.27), i.e. they are exact integrals of the generalized kinetic equations. 

Expressions (3.41) and (3.42) for 4-current end energy-momentum tensor, respectively, 
as well as self-consistency relations (3.47) and (3.48) are written explicitly for the case 
of nonrelativistic particles which number is conserved. This follows from the conventional 
way of nonrelativistic renormalization for such particles based on normal ordering. When 
the number of particles is not conserved (e.g., for phonons) or a system of relativistic 
particles is considered, one should replace F(X,p) — > | (^F(X,p) =f F(X,p)^j in all above 
formulas in order to take proper account of zero point vibrations (e.g., of phonons) or 
of the vacuum polarization in the relativistic case. These symmetrized equations, derived 
from special (=f) combinations of the transport equations (3.27) and (3.28), are generally 
ultra-violet divergent, and hence, have to be properly renormalized at the vacuum level. 
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3.5 Positive Definiteness of Kinetic Quantities 



For a semi-classical interpretation one likes to have the Wigner distributions F(X,p) and 
F(X,p) to be positive semi-definite, hereto after just called "positive". Using the operator 
definition for the Green's functions (2.2) and integrating it over a large space-time volume, 
one arrives at, e.g., 



J dXF(X,p) = ((J dye*V( y )) (/ dxe'^^x)^ , (3.49) 

and similarly for self-energies £ expressed through the current-current correlator. It in- 
dicates that the r.h.s. of such equality is real and nonnegative. Thus, we get the following 
set of constraints 



J dXF(X,p) > 0, J dXF(X,p) > 0, j dXT out (X,p) > 0, J dXT- m (X,p) > 0. (3.50) 

Similar relations are obtained for the integration over four-momentum space rather than 
space and time. As a result, in stationary and spatially homogeneous systems, in particular 
in equilibrium systems, the quantities F, F, T in and r out are real and nonnegative, i.e. 



F(p)>0, F(p)>0; r in (p)>0, r out (p)>0. (3.51) 

In deriving constraints (3.50) and (3.51), we did not use the fact that the Green's func- 
tions are solutions of the Dyson's equation. However, we used the operator picture. Any 
approximation, in particular, if formulated in the space of Green's functions, may spoil 
such rigorous statements like (3.49). Nevertheless, both the ^-derivable scheme and the 
gradient approximation preserve the retarded relations (A. 6) among the different con- 
tour components and the retarded and advanced functions of any contour function, with 
definite values for the imaginary parts of the corresponding retarded Wigner functions 



- 2lmG R (X,p) = A(X,p) > 0, -2lmE R (X,p) = T(X,p) > 0, (3.52) 

which even hold locally. In particular, solution (3.21) for the retarded Green's function 
shows that all retarded relations hold locally: the momentum part is the same as that 
in the homogeneous case with the space-time coordinate X as a parameter. Under the 
condition \ f — 7| <§; 1, cf. the discussion around eq. (3.24), one finds that 



r out (x, P ) » fT(x, P ) > 0, r in (x, P ) « (i T f)r(x, P ) > o (3.53) 

as long as the Wigner densities / and 1 =F / are positive. As the gradient approximation is 
a quasi-homogeneous approximation, one may therefore expect the positivity of T in and 
r out to be preserved even in the self-consistent treatment discussed here. Diagrammatic 
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rules may also corroborate this, since diagrams for r in and r out are calculated like in the 
homogeneous case. 

We now like to show that, if r in and r out are positive, also under minor restrictions 
the kinetic equations (3.27), (3.28) preserve the positivity of F and F, once one has 
started from positive F and F initially. For this purpose we turn to the generalized kinetic 
equation (3.27) for the phase-space distribution F(X,p) which is of integro-quasi-linear 
first-order partial differential type in 8 dimensions. We present it in the form 



PF = Ar in -rF+|r^,ReG fl |, (3.54) 

with the coefficient Z (cf. eq. (3.32)) in front of the time derivative being positive. The 
condition Z > takes place due to the dispersion relation for the self-energy S R and 
positive definiteness of F(X,p), cf. eq. (3.52). 

Given the solution, let us then discuss properties of F along characteristic curves deter- 
mined by the quasi-linear drift operator T> with a curve parameter s growing monotonically 
with time. Assuming F is positive initially, the occurrence of a negative value at some 
later time requires the first zero value to occur at one of the characteristic curves. Be s 
such a place, one finds 



s = s 



II V-.HvC 11 



s = s 



(3.55) 



If the r.h.s. is positive, one comes to a contradiction, as the approach from some positive 
value to a zero value which then becomes negative would require a nonpositive deriva- 
tive at s - Thus, AT in > \ {VF/A, ReG R j | is a sufficient (not necessary) condition that 

F preserves its positivity. The corresponding condition AT out > \ {VF/A, ReG^j | pro- 
vides positivity of F. According to eq. (3.30) this puts restrictions on the size of the 
occurring gradients of the gain and loss rates along the characteristic curves, i.e. along 
the classical paths. This condition is met in a wide range of cases, namely (a) in the 
quasi- homogeneous and quasi-static case, i.e. close to equilibrium, where gradients are of 
order C <C AT in , AT out , and (b) in the standard quasi-particle limit, where an extra small 
parameter, namely the width T, enters the estimations. 



4 Collision Term 



To further discuss the transport treatment we need an explicit form of the collision term 
(3.16), which is provided from the <£> functional in the — h matrix notation via the variation 
rules (A. 8) as 
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C(X,p) = -J- F(X,p) - — — , F(X,p). (4.1) 

Here we assumed $ be transformed into the Wigner representation and all =FiG~ + and 
iG H to be replaced by the Wigner-densities F and F . Thus, the structure of the collision 
term can be inferred from the structure of the diagrams contributing to the functional $. 
To this end, in close analogy to the consideration of ref. [17], we discuss various decomposi- 
tions of the ^-functional, from which the in- and out-rates are derived. This consideration 
is based on the standard real-time diagrammatic rules, where the contour integrations are 
decomposed into two branches with — and + vertices for the time- and anti-time-ordered 
branches, cf. Appendix A. For the sake of physical transparency, we confine our treatment 
to the local case, where in Wigner representation all the Green's functions are taken at 
the same space-time coordinate X and all nonlocalities, i.e. derivative corrections, are 
disregarded. Derivative corrections give rise to memory effects in the collision term, which 
will be analyzed separately for the specific case of triangle diagram (see subsect. 4.4). 



4-1 Diagrammatic Decomposition into Physical Sub-Processes 



Consider a given closed diagram of $, at this level specified by a certain number n x of 
vertices and a certain contraction pattern which links all vertices with lines of certain 
arrow sense for complex fields. This fixes the topology of such a contour diagram. It 
leads to 2™ A different diagrams in the — h notation from the summation over all — h signs 
attached to each vertex. Thus for any — h type diagram of $ also the diagram, where all + 
and — vertex signs are interchanged, contributes to $. Furthermore, for any diagram with 
given line senses the diagram with all line senses reversed also is a valid diagram of $. We 
shall exploit these two discrete symmetry operations to further determine the properties 
of <E>. The simultaneous application of the interchange of all vertex signs and the reversion 
of the line sense leads to the adjoint expressions, since in the underlying operator picture 
it adj ungates all operators and inverts the operator ordering. The corresponding values 
are then complex conjugate to one another, cf. (A. 7). 

Using these symmetries one can convert the functional $ into the following general form 



1$ = ^ ^ + i$ m',m) , ( 4 - 2 ) 

m,m' 

_rd*p ] _ dVd^ d^ni^fc 
iyW-y ( 27T )4 ' ' ' ( 2vr )4 (2vr) 4 " ' (2vr) 4 [^ Pt 

x R mM (X;p 1 ,...,p m ;p' 1 ,...,p' m ,)F 1 --- F m F[ ■ ■ ■ F' m , , (4.3) 

where in view of the local approximation the four momentum conservation has been 
extracted. While $ o compiles with terms void of any Wigner densities, i.e. from diagrams 
where all vertices have the same sign and which do not contribute to the collision term, the 



20 



nontrivial $ m , m ' terms sum the sub-class of diagrams of $ with precisely m + m' Wigner 
densities Fi = F(X,pi) and F[ = respectively. According to eq. (4.1) each $ m , m ' 

in (4.3) generates a multi-particle gain or loss contributions expressed in terms of integrals 
over products of generalized distribution functions F, and Fermi/Bose factors F. Every 
term in the sum (4.2) has been duplicated, repeating each term in its line reverse form. 
The corresponding transition rates R m ^ m i{X; pi, . . . ,p m ]p[, . . . ,p' m ,) are real due to the 
adjungation symmetry. As explained below, they result from the product of chronological 
Feynman amplitudes, given by the sub-diagram compiled from all — vertices linked by 
iG Green's functions times that of the anti-chronological part containing the iG ++ 
functions. 

The generic form (4.3) of $ m , m ' can be illustrated diagrammatically. For this purpose 
consider a — h notation diagram contributing to $ m>m ', which contains vertices of either 
sign. It can be decomposed into two pieces, denoted by a compact bra-ket notation, say 
(a | and \/3), in such a way that each of the two sub-pieces contains vertices of only one 
type of signf^] 



i$ 



a/3 



n, 



a/3 




— (a F 1 ---F m F[---F' m , 



— J (2^)4^0 6 (J> " 

x V*{X; p mm ,) F\ - ■■ F m F[ ■ ■ ■ F' m , Vp{X; p r , 



(4.4) 



where n a p counts the number of vertices in the closed diagram (a\/3). The short-hand 
notation p mm ' = {pi, . . . ,p m ;p[, . . . ,p' m >} summarizes the momenta, type and internal 
quantum numbers of the set of ordered valences, to be joint to the product of m + m! 
Wigner densities^ F\ ■ ■ ■ F m F[ ■ ■ ■ F^ linking the two amplitudes. The "end-caps" (a| and 
|/3) represent multi-point vertex functions, in simple cases of tree type, of only one sign 
for the vertices, i.e. they are either entirely time ordered (— vertices) or entirely anti-time 
ordered (+ vertices). Each such vertex function 



\P) = Vp(X;p mm ,), (a\ = V:(X;p mm/ ), (4.5) 

to be determined by normal Feynman diagram rules, has m F-valences and m' F-valences, 
respectively. In (4.5) we used the fact that adjoint expressions are complex conjugate to 

3 To construct the decomposition, just deform a given mixed-vertex diagram of $ in such a way 
that all + and — vertices are placed left and respectively right from a vertical division line and 
then cut along this line. 

4 This product of Wigner-densities originates from the =FiG ^ and iG^ Green's functions, 
(cf. eq. (3.7)). In closing the diagram by these Wigner densities extra fermion loops may appear 
besides the ones accounted for in the amplitudes a and j3. Using the product F\ ■ ■ ■ F m F[ ■ ■ ■ F' m , 
just take care of the additional (—1) factors arising from this closing of the diagram. 
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each other, cf. (A. 7). Accumulating all diagrams of $ that lead to the same set of Wigner 
densities F\ • • • F m F[ ■ ■ ■ F^ provides us with the generic form (4.3) with the partial process 
rates 

Rm,m'(X',Pmm') = E —Re {V*(X ; p mm ,)Vp(X ; p mm/ )} . (4.6) 

(°0)e* m w nQ/3 



The restriction to the real part arises, since with (a\j3) also the adjoint {(3\a) diagram 
contributes to the subclass $ m , m '. However these rates are not necessarily positive as in 
perturbation theory^. In this point, the generalized scheme differs from the conventional 
Boltzmann kinetics. Even the two rates, R mm > and its line sense reverse R m 'm, have not to 
be necessarily identical to each other. Still eq. (4.3) represents the most general form of $ 
expressed through the Wigner densities F(X,p) and F(X,p) in the local approximation. It 
is however important to realize that in many physically relevant cases, e.g. those discussed 
below, one indeed finds that 



Rm,m' Rm',m- (^"'O 

This property will be used as a sufficient condition for the derivation of the H-theorem. In 
the following we will restrict the discussion to cases where (4.7) is assumed. The treatment 
of more general cases, in particular in connection with the H-theorem, will be deferred to 
another publication. 

Since $ is two-line irreducible, there are at least three lines connecting (oj| and \(3) and in 
many cases (a| and \(3) are connected by, at least, four lines, like in the case of two-body 
interactions. In such diagrams each of the amplitudes a or j3 necessarily form a connected 
diagram for the complex field case for binary (m = m! = 2) and for triple scattering 
(m = m' = 3) case, while in the general case of multi-particle (more than triple) processes 
such amplitudes may be disconnected. 

No further symmetry can be specified at this level without additional knowledge on possi- 
ble topological and other symmetries of $ and the different particle species involved. The 
decomposition discussed here solely relies on a straightforward application of the contour 
rules for multi-point functions. They differ from other cutting rules like those derived by 
Danielewicz [29], which represent the result in terms of generalized retarded functions. 



5 In perturbation theory they are positive, since there the sum in (4.6) leads to absolute squares. 
In the general case with resummed propagators this positivity may be lost due to the restriction 
to diagrams of which are globally two-particle irreducible. The latter excludes certain combi- 
nations of amplitude diagrams, implying that the rates of genuine multi-particle processes are 
not necessarily positive. 
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4-2 Local Collision Term and Memory Corrections 



The gradient corrections to the folding of the self-energies with the propagators in the 
collision term have already been accounted for in the kinetic equations. They essentially 
give rise to in-medium corrections of the convective part (l.h.s) of the generalized kinetic 
equation (3.27). Within the spirit of the gradient approximation one also likes to express 
the self-energies themselves by space-time local quantities. Thus, for a consistent gradient 
approximation further gradient corrections are admissible once the self-energy diagrams 
and thereby the diagrams of $ consist of more than 3 vertices. We call the collision term, 
evaluated with all Green's functions in the Wigner representation taken at the same 
space-time point X, the local collision term. 

In terms of the representation (4.3) of $ and implying line-sense reversal symmetry for 
the rates (4.7) the matrix variation rules (4.1) determine the following local part of the 
collision term (3.16) for a particle of flavor "a" as 



a [ ,P) ~ 2 A, J Wf ' ' ' (2?r) 4 (27T) 4 ' ' ' (2VF K ™,m>{*,Pl, • • • ,Pm,Pi, • • ' ,P m >) 

{Fi ■ ■ ■ F m F[ ■ ■ ■ F' m , - Fi ■ ■ ■ F m F[ ■ ■ ■ F' m , } 



m,m' 



X 



X 



S a ai S\pi -P)~Y1 6 a ai S\p'i - P) 
i=l i=l 



i=i 



The functional variations of $ with respect to F a and F a are expressed in terms of the 
four- momentum projectors S aai 5 4 (pi — p) which in case of multi-component systems of 
particles with different flavors and internal quantum numbers "a" also include the proper 
projections onto the different flavors = a. This expression nicely visualizes the detailed 
balance property namely that the same multi-particle rate determines both the forward 
and the backward processes. 

The first-order gradient corrections to the local collision term (4.8) we call memory cor- 
rections. We refrain from specifying the latter in general terms. They are linear in the 
spatial gradients of Wigner functions and the displacement factors in front of these Tay- 
lor expansions give momentum gradients of associated other multi-point functions. An 
explicit example will be discussed in section 4.4. Often (but not in general) there are 
special reasons why C mem is small. In particular, only diagrams of third and higher order 
in the number of vertices give rise to memory effects, which are discussed below (sect. 4.4) 
within a simple model. One has to keep in mind that memory effects are quite commonly 
neglected in transport models. 



23 



4-3 $ -Derivable Collision Terms 



4-3.1 Two-body Potential Interaction 

To be specific we consider a system of fermions interacting via a two-body potential 
V = V 5(x — y), and, for the sake of simplicity disregard its spin structure, by relating 
spin and anti-symmetrization effects to a degeneracy factor d. To derive the decomposition 
of a ^-derivable collision term, we employ the rules as described in subsect. 4.4. 

In the first example, we consider the generating functional $ to be approximated by the 
following two diagrams 




with the dashed line illustrating the decomposition according to (4.4). Here the l/n\ 
factors start with 1/2, 1/4, . . . according to the nonrelativistic diagram rules for two-body 
interaction, cf. Appendix C, i.e. the vertex dots are considered as the zero range limit 
of a finite range interaction. In the { — h} matrix notation of the Green's functions, one 
can easily see that the one-point diagram does not contribute to the collision term, while 
decomposing the second diagram along the dashed line leads to a purely local result 



r (2) = r2 f d V d 4 p 3 

J (2tt) 4 (2tt) 4 (2vr) 4 

x 5 A (p + Pl - p 2 - p 3 ) (F 2 F 3 FF 1 - F 2 F 3 FF 1 ) , (4.10) 

where the brief notation of the previous subsect. is used for Fj and Fj. This collision inte- 
gral has precisely the form of the binary collision term of Boltzmann-Uehling-Uhlenbeck 
(BUU), except for the fact that distribution functions are not constrained to the mass 
shell. The binary transition rate 




4 2) = v? = | x f ( 4 - n ) 

is nonnegative in this case. Here and below, the bracketed superscript ^ (or < - 3 - ) ) points 
out the origin of the quantity (C, R, etc.) from the second (or third) diagram of <3>. The 
subscript 2 in the transition rate of eq. (4.11) indicates the binary-collision nature of this 
transition rate. Note that external 4-momenta (in-going and out-going) of the scattering 
amplitude X are n °t confined to the mass shell. For the trivial case under consideration, 
this fact does not give rise to any important consequences. However, for more complicated 
examples below, it means that the collision term is determined by off-shell scattering 
amplitudes. 
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The picture becomes more complicated, if $ involves diagrams of higher orders. For in- 
stance, let us add the following three point diagram to $, i.e. 



i$ = i($(l) + $(2) + $(3)) 




(4.12) 



where one possible decomposition is illustrated by dashed lines. The corresponding self- 
energy becomes 



iS = -i (£« + £^ + £< 3 )) = 




(4.13) 



Now the collision term contains a nonlocal part due to the last diagram. This nonlocal 
contribution is discussed in the next subsect. in detail. The local part can easily be derived 
in the form 



+ cg> = d 2 



d 4 pi d 4 p 2 d 4 p 3 
(2tt) 4 (2tt) 4 (2tt) 4 



X + 



INI 



x 5 4 (p + Pl - p 2 - p 3 ) (F 2 F 3 FFi - F 2 F 3 FF l ) , 



(4.14) 



where all the vertices in the off-shell scattering amplitudes are of the same sign, say " — " 
for defmiteness, i.e. there are no " H — " and " — h" Green's functions left. The quantity 
(7(2) _|_ is again of the Boltzmann form with the transition rate 



„(2) „(3) 
K 2 + ri 2 



X + 



(4.15) 



At the first glance, one may argue that this rate is not necessarily positive in the limit of 
strong coupling. Indeed, the first term in eq. (4.15), i.e. — iVo, is purely imaginary, whereas 
the second one — the loop — given by 



loop = —d J 



-j—^ | V I 2 iG— {x,p + px)\G~~ 
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has both real and imaginary parts. Hence, the real part of this loop is canceled out in eq. 
(4.15). If | Im(loop) |>| V I, the rate {R> 2 ! + R K 2 >) may become negative, depending on 
signs of Vq and Im(loop). However, one has to keep in mind that the Green's functions in 
the loop cannot be chosen arbitrarily. Rather in a consistent treatment, as shown in eq. 
(80) of ref. [17], the loop reveals a factor oc |l/r| oc |1/Vq 2 | resulting from the imaginary 
part of the retarded self-energy in the propagators, which balances the total value of 
the loop. Indeed, in equilibrium the gain part of collision term and thus (R 2 + R 2 ) is 
positive. This illustrates that the question of positive definiteness is quite subtle. 



4-3.2 Bosonization of the Interaction 

It is obvious that the situation becomes complicated once one extends the picture to 
ring diagrams with more than three vertices. Yet, there is a simple (however, not gen- 
eral) strategy to proceed. We may avoid the two-particle potential interaction from the 
very beginning and rather introduce an interaction mediated by an artificially introduced 
neutral heavy scalar boson (h.b.) of mass mh.b. much larger than any characteristic mo- 
mentum transfers in the system. Then the free retarded Green's function of the boson 
approximately equals 



A 0R 



-1 



m 



h.b. 



-iO 



(4.16) 



and the vertex of fermion-heavy-boson interaction being g = y |Vo|mh.b.- Moreover, by the 
same reason the heavy-boson occupation numbers may be put to zero, = 0, since 

typical excitation energies are much less than the boson mass. The fact that this boson 
is very heavy makes the fermion-fermion interaction almost point-like. To be specific, we 
assume that Vo < 0, i.e. attractive, which can be mediated by a scalar boson. In case of 
repulsive interactions, a vector boson would be an appropriate choice. 

Thus, from now on, we deal with a system of interacting fermions and heavy bosons. Let 
us take the following approximation for the corresponding $ functional (we call it $h.b.) 




(4.17) 



in terms of full Green's functions of fermions (the bold solid lines) and bosons (the bold 
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wavy line). In this approximation, the boson self-energy is given by (cf. eq. (2.9))0 



.1„ <5i$h.b. 

2 <3lGh.b. 



(4.18) 



and the heavy-boson Green's function is defined by the standard Dyson's equation 



iCh.b. 



+ 




(4.19) 



We would like to compare this model with the model described in subsection 4.3.1. Elim- 
inating the artificial heavy boson one effectively sums up all ring diagrams of the type 



i$ 



ring 




(4.20) 



Note that $ r i ng ^ $h.b., since the summations of loops in eq. (4.17) and eq. (4.20) have 
different sense. In eq. (4.19) we have the conventional diagrammatic summation, while in 
eq. (4.20) the summation is logarithmic, i.e. with the factors 1/n, where n is a number 
of vertices in a diagram. However, this difference in summations is compensated for by 
the heavy-boson contribution to the generating functional (2.8). The equivalence of ap- 
proximations (4.20) and (4.17) can be actually seen, e.g., from the fact that they result 
in the same approximation for the fermion self-energy iEf = 5^} T - mg /5G{ = SQ^J SGf (cf. 
eq. (2.9)) after substituting eq. (4.19) for the heavy-boson propagator. Here and below, 
the sub- label "f denotes fermion quantities. 

Upon introduction of the heavy-boson degree of freedom, we have to deal with two cou- 
pled transport equations — for nonrelativistic fermions and for heavy bosons — with the 
following collision terms 



C ^ P)= J (2^(2^ 96 (P-^"») 



(2tt) 4 (2tt) 4; 

F i {X, P )F i {X, Pl )F h , h ,(x,p 2 ) ~ F i (X,p)F { (X, Pl )F h . h Xx,p 2 ) 



(4.21) 



for fermions and 



6 Note that for neutral bosons, which number is not conserved, the additional factor 1/2 appears 
in eq. (2.9), cf. ref. [8]. 
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C h . h {X.p) = d 




x 'F h . h .(X,p)F { (X,p 1 )F i (X,p 2 ) - F h . h .(X,p)F i (X,p 1 )F i (X,p 2 )] , (4.22) 



for heavy bosons, where g 2 = — Vqw^ b > is defined in terms of two-particle interaction 
strength Vq < and the heavy boson mass mh.b.- 

The collision terms (4.21) and (4.22) are very simple in spite of the fact that they involve 
the whole series of ring diagrams. The corresponding gain and loss terms are positive and 
contain no memory effects, as they are hidden in the boson, while in the pure fermionic case 
already the triangle diagram gives rise to memory effects. Indeed, there is no contradiction 
here. If one wants to eliminate the bosonic degree of freedom, one has to resolve the bosonic 
transport equation with respect to the bosonic generalized distribution function F^ for 
entire past and substitute this into the fermionic collision term. In this way, the resulting 
collision term becomes highly complicated and nonlocal and thus contains memory effects. 

Hence, we have demonstrated that sometimes it is useful to introduce new degrees of 
freedom in order to achieve a reasonable collision term. Of course, all the consideration 
above remains valid also for particle-particle interaction mediated by a real boson rather 
than only by artificially introduced one. 

4-4 Memory Effects in Collision Term 

A general treatment of memory effects in the collision term is a cumbersome task. In this 
subsection we continue to consider a system of nonrelativistic fermions interacting via the 
contact two-body potential and concentrate on the third diagram in eq. (4.12), which has 
already been considered in a ^-derivable scheme for the thermodynamic potential and 
entropy in refs [14,15]. This is the first ring diagram to contribute to memory effects. The 
corresponding self-energy diagram reads 



Standard diagrammatic rules in the matrix representation present E^- in analytic form 

as 



z 




(4.23) 



— i 



-(3) 
•jk 



dz ■ iV a jj >a kk riG : >' k ' 



(x,y)L k ' l (y,z)a ll/ L l ' j ' \z,x), 



(4.24) 



where we have introduced a loop function 
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L jk (x, y) = fe <*~> = di mCP k (x, y)iG*(y, x), (4.25) 

and Oij are given by eq. (A. 2). As above, the factor d results from the trace over spin. In 
the Wigner representation, U k takes the form 

# W) = / ^V\X-y+j,J), (4.26) 
where 

P k (X;p"+p',p") = diV iG jk (X,p"+p')iG kj (X,p"). (4.27) 
The loop functions L possess notable properties 

L ++ + L-=L-+ + L+-, Li k (x,y) = L k i(y,x), L jk (X, p) = L kj (X, —p) . (4.28) 

The former property follows from the general property (A. 6) of the two-point functions 
and holds in both the coordinate and Wigner representations. Proceeding from eq. (4.24) 
and with the help of relation (3.5), we can immediately evaluate the Wigner transform of 
(cf. eq. (3.1)) and perform its gradient expansion 

Eg. (E-L +(£<?)_, (429) 
where 

-i(Eg ) ) ioc (X,p) = J |^iy ^ fefe ,a^^ 
is the local contribution to S-^, and 



- 1 (^L. ^ = \ I ^^ f a kk/ a ul iG^(X,p'+p) 



X 
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{L k ' l (X,p'),L l 'i'(X,p')}^ x (4.31) 



is the first-order gradient (memory) correction, which is of central interest in this subsec- 
tion. Here, the Poisson bracket is taken with respect to the (jpf , X) variables indicated in 
the subscript. 

To evaluate the memory correction for the collision term (cf. eqs (3.9) and (3.16)), we 
need to consider the following combinations of self-energies with Green's functions 
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( E - )_ (*,P)G-+(X,p) = i / {L-', !.'+} ,(4.32) 



(2ir) 4 (i 



<7-(x,rt (s™ )_ (A',p) = 1 / ^ir*d;p' + P,p)* {i+<. L'-}^ . (4.33) 
With the symmetry relations 

= <*<* = {^.^U (4 - 34) 

deduced from (4.28) one determines the first-order gradient correction to the collision 
term induced by graph (4.23) as 



C$L(X,p) = [(S^) mem (X, P )G-+(X, P ) - G+-(X, P ) (E?l) mem (X, P ); 



i /■ d 4 p' 1 
~ 2 J (2vr) 4 rf 

The corresponding local collision term is given by eq. (4.14). 



5 H— Theorem 



5. 1 Time-Irreversibility of the Generalized Kinetic Description 



Compared to exact descriptions, which are time reversible, reduced description schemes 
in terms of relevant degrees of freedom have access only to some limited information and 
thus normally lead to irreversibility. Various reduction and coarse graining schemes have 
been discussed and developed over the years. In the Green's function formalism presented 
here the information loss arises from the truncation of the exact Martin-Schwinger hier- 
archy, where the exact one-particle Green's function couples to the two-particle Green's 
functions, cf. [6,22], which in turn are coupled to the three-particle level, etc. Q This 
truncation is achieved by the standard Wick decomposition, where all observables are 
expressed through one-particle propagators and therefore higher-order correlations are 
dropped. This step provides the Dyson equation (2.3) and the corresponding loss of in- 
formation is expected to lead to a growth of entropy with time. 

At the operator level the nonequilibrium entropy can be straightforwardly formulated in 
terms of von Neumann's entropy, which is an entropy in the information theory sense [30] 



7 The appropriate classical limit leads to the Bogolyubov-Born-Green-Kirkwood-Yvon chain 
of equations. 



30 



S = -Trpln p. 



(5.1) 



It is given by the expectation value of the logarithm of the density operator p itself. The 
problem of the Green's function formalism is that it does not give a direct access to the 
density operator itself but rather describes the space-time dependence of the expectation 
values of well defined operators, say (^A(t)^ = Tr A p(t). To this extend, there is no im- 
mediate formulation of the entropy (5.1) in terms of Green's functions for nonequilibrium 
case. The situation is different at equilibrium, where with p = exp(—j3(H — pN))/Z a 
well defined density operator exists, which leads to the well known Matsubara or real- 
time formulation of the equilibrium entropy, cf. sect. 6 below. As will be seen, even in 
equilibrium, the entropy expression is not priori given but rather depends on the choice 
of $ in a ^-derivable scheme. 

In order to access a nonequilibrium expression relevant for our generalized transport 
equation (3.27), we shall start from this transport equation and derive a flow expression 
s M with the property that its divergence grows in time, i.e. 



> 0, (5.2) 

and which in the equilibrium limit merges the corresponding equilibrium form of the 
entropy flow. That is, we show the existence of an f/-theorem for our generalized kinetic 
description and s^(X) is thus identified with the kinetic entropy flow. Thereby, the validity 
conditions for the derivation of this kinetic entropy flow coincide with those of the kinetic 
equations themselves. 



5.2 Markovian Entropy Flow 



We start with general manipulations which lead us to definition of the kinetic entropy 
flow. We multiply eq. (3.28) by (=f) \n(F/A), eq. (3.27) by -\n(F/A), take their sum, 
integrate it over d 4 p/ (27r) 4 and finally sum the result over internal degrees of freedom like 
spin (Tr). Using the identity for the Poisson brackets 

{B, Af} In / ± {B, A(l T /)} In(l T f) = {B, Af\nf± A(l T f) m(l =F /)} (5.3) 
for any functions A, B and positive / and 1 =F /, one then arrives at the following relation 



where the quantity 
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>loc 



'loc, a 



TrE 



d 4 p 



<9Re£^\ 



ReGf Tin 



Fa_8_ 

A a <% 





A a dp^ V a A n 




(5.5) 



obtained from the l.h.s. of the kinetic equation is interpreted as the local (Markovian) 
part of the entropy flow. Here we have restored the summation over "a" denoting the 
different particle species and intrinsic quantum numbers for a multi-component system. 
It illustrates that, although this entropy expression accounts for interactions among all 
particles, it can be expressed as a sum of the individual contributions, each of which is 
solely determined by the particle self-energy Re£„ and its width T a . Gradient corrections 
to the collision term (i.e. C mcm ) give rise to extra memory contributions to the entropy 
flow. 

Partial integrations in eq. (5.5) lead us to a more transparent expression for the entropy 
flow 



^Ej-^ A L(X,py a (X,p), (7 a (X,p) = =F[l=F/]Ml=F/]-/ln/, (5.6) 



'loc 



(2k) 



where 



(5.7) 



has the meaning of an entropy-flow spectral function, while 



— A 



P - 



<9Re£ 
dp. 



M a T- 



1 d£a 
dp^ 



(5.8) 



is the flow spectral function, cf. the corresponding drift term (proportional to d^f in 
eq. (3.24)). The zero-components of these functions coincide with the corresponding two 
functions introduced in refs [31,32] for the case of equilibrium systems. Moreover, the B° 
function satisfies the sum rule 



/ 



d Po R o 



(5.9) 



which is easily obtained from the sum rule (3.23) for the spectral function A with the help 
of Kramers-Kronig dispersion relations. For the case that the considered particle a is a 
resonance, like the A or p-meson resonances in hadron physics, the B° function relates to 
the energy variations of scattering phase shift of the scattering channel coupling to the 
resonance in the virial limit, for details see e.g. refs [31,32,8]. 
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In the noninteracting-particle limit, the entropy (the zero component of s{^ c ) directly 
transforms into the proper ideal gas expression, cf. ref. [33]. In the quasiparticle approxi- 
mation, the T in and r out terms have an additional smallness, which allows to neglect these 
terms. Thus, expression (5.5) for the entropy flow takes the form 



Koc) 



SO 



TrE 



f d 3 p 


f 1 ^ 




1 (2tt) 3 


de a 




V dp) 





T(lT/7)ln(l T / a qP )-/ a qP ln^ 



qp 



(5.10) 



in the quasiparticle limit, which follows from the substitution of eqs (3.35) and (3.36) 
into eq. (5.5). From (5.6) eq. (5.10) is also easily recovered, since A 2 T /2 transforms to the 
corresponding 5-function in the limit r — > 0. 

To prove explicitly the if-theorem we have to show that the r.h.s. of eq. (5.4) is nonnega- 
tive. To this end, we should consider the convolution of the collision term with \n(F/F). 
First, we do this for collision terms in local approximation. 



5.3 Local Collision Term and H-Theorem 



Using the multi-particle process decomposition (4.8) we arrive at the relation 

d 4 pi d 4 Pm d 4 p[ d A p' m , 



Tr/ 7 ^ln^ 1 oc(X,p) = TrE 1 /^ 

J (27f) 4 p loci ,P) A. 2 J 



X 



(2vr) 4 F ^ ^2 J (2tt)' (2n)*(2«)* (2k) 

{F^■■■F F[---F' ,-Ft--F F[ ■ ■ ■ F' ,\ In ' ' ' 1 ' ' ' m ' 



4 



v i=l j=l 



Here we assumed different flavors and intrinsic quantum number to be absorbed in the 
momenta p\ and p\. In the case when all rates R m ,m' are nonnegative, i.e. R m , m ' > 0, this 
expression is nonnegative, since (x — y) ln(x/y) > for any positive x and y. In particular, 
Rm,m> > takes place for all $-functionals up to two vertices. Then the divergence of s^ oc 
is nonnegative which proves the if -theorem in this case with (5.5) as the nonequilibrium 
entropy flow. 



5.4 Explicit examples for the H-Theorem 



We explicitly discuss the two examples introduced already in sect. 4.3 with Markovian 
collision terms, i.e. with the $ functional consisting of one- and two-point diagrams only. 
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In the pure fermionic case with collision term (4.9), one can state an exact if -theorem 

r d 4 p F {2) 3 r d 4 p d 4 p x d 4 p 2 d 4 p 3 i 12 
J (2tt) 4 F J (2tt) 4 (2tt) 4 (2tt) 4 (2tt) 4 I I 

.P.?" 1 F F ~ ~ ~ ~ 

x5 4 (p + Pl -p 2 - p 3 ) In (FF1F2F3 - FFiF 2 F 3 ) > 0. (5.12) 

FF 1 F 2 F 3 v 7 

Note that here and below the extra factor 0? comes from summation (Tr) over spin of 
fermions. 

Furthermore, it is instructive to consider the $ approximation (4.17), where a heavy scalar 
boson has been introduced in order to sum up the series of ring diagrams of eq. (4.20), cf. 
subsect. 4.3.2. In this case, there are two coupled kinetic equations with collision terms 
(4.21) for fermions and (4.22) for heavy bosons, leading to 

«c = 0X + «.b. = / T^fl jC* + 5 ^ §7^b.) , (5.13) 

where s{^ c is given by the sum of the proper fermion (s^ ) and heavy-boson (s^ ) contri- 
butions. The thermodynamic entropy for this system, i.e. for $ given by diagram (4.17), 
has recently been obtained in ref . [39] . Here we present the corresponding nonequilibrium 
entropy flow together with an affirmation of the if-theorem also for this case. The collision 
term, indeed, becomes 

o fj, , f d 4 p 1 d 4 p 2 d 4 p b 2 4 

d ^ = d J (2^(2^(2^^ (Pi-P2-Po) 

x In l lF ~ F ~ h - (FflFfiFh.b. - F n F f2 F h . h ) > 0, (5.14) 

ffl^ f2^h.b. 

which is nonnegative. 

Our representation of the entropy of a system interacting via two-body potential (i.e. as a 
sum of a purely fermionic part and that of the artificially introduced heavy boson) is also 
very similar to that derived by Riedel [14] within the ring-diagram model of ^-derivable 
thermodynamics. In both cases, the bosonic part of the entropy s^ h takes account of the 
fermion-fermion interaction calculated within the ring-diagram approximation (4.20). In 
thermodynamics, this interaction part of the entropy gives rise to the famous correction 
to the specific heat of liquid 3 He [14-16]: ~ T 3 lnT, where T is the temperature. As has 
been found by Carneiro and Pethick [15], this correction to the specific heat emerges 
already solely from the third diagram of the whole ring series (4.20). To demonstrate 
the same within our kinetic approach, we should consider the ^-derivable model (4.12) 
involving only the first three ring diagrams. Moreover, since the local entropy expression 
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(cf. eq. (5.5)) derived above does not contain such kind of corrections, one has to explicitly 
consider memory effects in the collision term (4.35). 



5.5 Memory Effects in Entropy Flow and H-Theorem 

We assume that the fermion-fermion potential interaction is such that the corresponding 
transition rate (4.15) is always nonnegative, so that the if -theorem takes place in the 
local approximation, i.e. when we keep only + Cj^.. Our aim now is to derive the 
entropy, which takes into account memory effects in the collision term (C^ m ). 

Proceeding similarly to that in subsect. 5.2, we multiply eq. (3.28) by — \n(F/A), eq. 
(3.27) by —hi(F/A), make their sum and integrate it over d 4 p/(2-7r) 4 . Thus, we arrive at 
the equation 

McPO = Tr / yj^ In j{C^ + C£) + Tr J ^ In ^c£L, (5.15) 

where s^ oc is still the Markovian entropy flow defined by eq. (5.5). Our aim here is to 
present the last term on the r.h.s. of eq. (5.15) in the form of full x-derivative 

Tr J ^ In ^Cgl m = -d^ mcm (X) + 5c mcm (X) (5.16) 

of some function s mcm (X), which we then interpret as a non-Markovian correction to the 
entropy flow of eq. (5.5), plus a correction (5c mem ) which is small in some sense. Indeed, 
this term on the r.h.s. of eq. (5.16) is linear in X- and p-derivatives. Hence, it cannot be 
transformed into sign-definite form. The only possibility remained is to construct a full 
derivative of it. If we succeed to find a proper s mem (JT), then relying on smallness of 5c mcm 
we obtain 

d, « c + * me J * Tr J ^ In ^(C^ + C© > 0, (5.17) 

which is the if-theorem for the non-Markovian kinetic equation under consideration with 
s ioc + s mem as the proper entropy flow. The r.h.s. of eq. (5.17) is nonnegative due to our 
assumption that the corresponding transition rate (4.15) is always nonnegative. 

Hence, let us consider the last term on the r.h.s. of eq. (5.15). Upon substituting expression 
(4.35) for C^2 m and shifting the integration variable p — > p — p'/2 in this term, we arrive 
at 
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Tr 



(27T) 



L + -(X;p+ £,p - ^) - Z-+(X;p+ ^,p - 



x 



2 J (27r) 4 (27r) 4 

| L+- (X, p') , L~+ (X, p') } In * * 
I v v ,FJ i P ',x F(X,p-p> 



F(X,p-p'/2) 
/2) 



(5.18) 



By making the half-sum of eq. (5.18) with that with replaced p' —p', we arrive at the 
symmetric form of this equation 



d 4 p 

W) 



Tr/^ln^ 3 ) =1 

7 (2tt) 4 F mem 4 7 (2tt 



i f d 4 p d 4 p' 



(2tt) 4 (2ttY 



Z~+(X;p+ |,p - |) - Z + -(X;p+ |,p - | 



x{L- + (X,p / ),^ + -(^,p / )V,xln 



L + -(X;p + p'/2,p-p'/2) 



L-+(X;p + p>/2,p-p>/2y 
which after simple algebraic transformations can be decomposed into two terms 



Tr 



d 4 p 



ln-C® 



mem ''mem T wo mem , 



(5.19) 



(5.20) 



with 



.P0 



d 4 p d 4 p' 



L-+(X,p'),L + -(X,p')L- 



x 



In 



(27r) 4 (2vr) 4 k 
Z + -(X;p + p72 !P -p72) 
L-+(X;p + p'/2,p-p'/2) 



X;p+ r v ~ 2 



(5.21) 



<5c mcm (X) 



i /" d 4 p d 4 p' 



2 J (2vr) 4 (27r) z 



(X, p>) , L + - (X, p') } p , x Z-+fX;p + |,p- 



p 



L- + (X,pQ,Z- + fx;p + ^,p-^ln Z+ "^^ + ^^-^ 2 n 
. 1 ^ V 2' P 2j L-+(X;p + pV2,p-pV2)/ p , 



x L + -(X,p') 



(5.22) 



Using partial integration, we have subdivided the quantity of eq. (5.20) in such a way 
that the first term c mcm takes the form of the full divergence, and thus defines the non- 
Markovian contribution to the entropy flow (cf. eq. (5.16)) 



= - if ^ Z-+ X;p + P - ,p - P - ) L+-(X,p') 

lV ; 2 J 2n 4 2tt 4 I ,F 2' F 2 K ,y ' 
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x ^ L+~(X;p + p>/2,p-p>/2) 
X * L-+(X;p + p>/2,p-p>/2) 

This c mem -term remains nonzero even in the local thermal equilibrium, which is defined 
by the former equilibrium relations (B.1)-(B.4) but with temperature T(X), 4- velocity 
U U (X) and chemical potential fi(X) depending on the coordinates. On the contrary, as 
we show below, the second term <5c mem vanishes in the limit of local thermal equilibrium. 



dL-+(X,p>) 
dp' u 



(5.23) 



In local thermal equilibrium the Kubo-Martin-Schwinger condition (B.l) provides the 
following relations 

f L+~(X;p + pi/2,p-p>/2) \ [ L+-(X;p>) \ _ ( p' v U»{X) \ 

U-+(*;P+P72,P-PV2)L eq . U- + (^p')J loc . eq . exp l nx) 



which can also be derived from (B.2)-(B.5) proceeding from definitions of L* J (4.26) and 
jjj (4.27). Guided by (5.24) we write the L u ratio in the ln-term of (5.22) as 

L+-(X-p + p'/2,p-p'/2) _ L+-(X;pf) 

L-+(X;p + p'/2,p-p'/2) L-+(X,p>) { U ' K ' > 



where £ is expected to be small within the validity range of the generalized kinetic equa- 
tion (3.27), i.e. |£| ~ |/ — 7]. Substituting this into expression (5.22) for <5c mcm , the 
ln(L~ + /L~ + )-term in the second Poisson bracket cancels against the first term, since the 
p-integration converts the linear L lj factor into L lj factor. Thus one obtains 



Sc mcm = -- [ -^--^4fL- + ,L- + ln(l + 0} L + - 
2 J (2tt) 4 (2tt) 4 I ' v W V,* 

Ji f d 4 p dV r + ~ 1 1 

~2J (2vr) 4 (2vr)U ' 2^ ) p ,, x ' 

where also the term linear in £ exactly cancels out, since 



i r d 4 p dV r + ~ + + _ 
~2 7(2^(2^1 ' 

As above, the L* J can be readily integrated over p to produce L % \ revealing a cancellation of 
the terms in the Poisson bracket in (5.27). Thus, <5c mcm is not only zero in local equilibrium 
(£ = 0). It is of the second-order in the small parameter |£| m \ f — 7I times gradients and 
therefore negligible compared to c mem in (5.21) which is of linear order in |£| \ f— j\ times 
gradients. Thus within the validity of the generalized kinetic equation (3.27) s^ em (X) 
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as given in (5.23) represents the appropriate non-Markovian memory correction to the 
entropy flow. In local equilibrium expression this term can be further simplified to 



'mem)eq 



lnQM-1 

L-+(X,p>) 



dL' + (X,p') 



(5.28) 



by means of relation (5.24). 



6 Thermodynamic Limit of Entropy 



6.1 Thermodynamic Entropy 



In the Matsubara technique, the thermodynamic potential Q (see, e.g., ref. [4]) is a func- 
tional of Matsubara Green's functions G(e n ,p). Its important property is that it is sta- 
tionary under variations of G(e n ,p) at fixed free Matsubara Green's function G°(e n ,p) 



5Q 



5G(s n ,p)J G0 



0. 



(6.1) 



Since G(e n ,p) has the spectral representation 



G(e n ,p) 



de' A(e',p) 
2n e n — e' ' 



(6.2) 



property (6.1) implies that fl is stationary under variations of the spectral density A(e,p), 
provided the Matsubara frequencies e n , over which the summation runs, are not changed 
and, thus, the free Matsubara Green's function G°(e n ,p) is unaltered. 

In the standard way, cf. (B.8)-(B.10), the e„-sum is converted into an energy integral 
over the distribution functions n(e), cf. (B.4), and thermodynamic potential expressed in 
terms of the real-time Green's functions and self-energies becomes [15] (cf. eqs. (B.2) and 
(B.3)), 



n(s) (-21m In [-G R (p + i0,p)] - ReG R T - AReS^) + $ T (6.3) 



where e = po — \x in the rest frame of the equilibrated system, and $t is represented by 
the same set of closed diagrams as $. The thermodynamic entropy density is defined as 



s = -(dn/dT) fi /v, 



(6.4) 
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i.e. as the temperature derivative at the fixed chemical potential /i with V being the 
system's volume. Due to the stationarity property of Q, the only T-dependences to be 
taken into account in calculating the entropy are the explicit dependence of the discrete 
Matsubara frequencies on T or in the occupations n(e) in the integral formulation (6.3). 
From eqs (6.3) and (6.4) one finds 



- Tr / 



d 4 p dn(e) 
(2tt) 4 dT 



-21m In 



-G R ( Po + iO,p) 



- ReG R T - AReS 



R 



d<5> 7 



(6.5) 



for the entropy density. With the help of the identity (B.6) we replace dn(e)/dT in eq. 
(6.5) and then perform partial integration over p . Thus, we obtain 



S = S } 



loc S me m; 

4„ 



where 
d 



-21m ]n[-G R (po + iO,p)] - ReG H T 



Ri 



<7 

'—'mom 



9$ 



d 4 p g(AReS fi ) 



+ Tr 



/ 



<7(e) = =f(1 =F n(e)) ln(l =p n(e)) - n(e) \nn(e). 



(6.6) 
(6.7) 

(6.8) 
(6.9) 



We have used subscripts "local" and "memory" to denote these two different contributions, 
because below we demonstrate that they are indeed associated with local (Markovian) Sj° oc , 
cf. eq. (5.5), and memory (non-Markovian) sj, em , cf. eq. (5.23), parts of the kinetic entropy. 
Taking derivatives in eq. (6.7), we readily get 



'loc 



Tr 



d 4 p 



<j(e)A° s (p) 



(6.10) 



with A® defined in eq. (5.7). Thus, S'loc, indeed, coincides with 0-component of the kinetic 
entropy flow (5.6). 

In order to clarify the meaning of the values S'loc and S mem we first inspect the quasiparticle 
limit, in which the spectral function is reduced to the delta-function, cf. eq. (3.33). In this 
limit, the value (6.10) for S\ oc is reduced to 

= Tr / ( £{P) ~ ") ' (6 - U) 

where the quasiparticle energy e(p) is determined by solution of dispersion equation (3.34). 
The full Sfo c is just the sum of single-particle contributions, as if we deal with a noninter- 
acting ideal gas of quasiparticles. This is a standard picture in the quasiparticle approxi- 
mation. Corrections to Sf^, resulting from Si oc , are of higher order in the width T. At the 
same time, S mem provides corrections to S^ c even in the zero-order in T, which are associ- 
ated with real rescatterings of on-mass-shell quasiparticles. This fact was demonstrated by 
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Carneiro and Pethick in ref. [15]. For the Fermi liquid, they showed that the first diagram 
for <3> that contributes to S mcm is the triangle one $ ( - 3 * ) of eq. (4.12). Contributions of 
and $( 2 ) are zero. In the quasiparticle approximation, the contribution of $^ 3 ^ is (cf. eqs 
(45) and (72) of ref. [15]) 



oqp _ _ 
mcm 12 J (2tt) 4 



d4 P (r^ p (p)) 3 drih ^ 



dT 



(6.12) 



where rib is the thermal occupation number of artificial boson, given by eq. (B.4), and the 
bosonic width T h is defined as 



ri{ \e' uj ) rif ( e' H — oo 



A { (p'-^A { (p'+^ (6.13) 



with e' = p' Q — \i and uj = p , n f stands for thermal fermion occupations, cf. eq. (B.4). 
As above, V stands for the strength of the two-body potential. To get the quasiparticle 
approximation to this width (r^ p ), we should replace the exact spectral function A{ in 
eq. (6.13) by its quasiparticle approximation A^ p defined by eq. (3.33). Note that S°® em 
is expressed in terms of bosonic quantities rib and T b , although we have started initially 
with the purely fermionic system with potential interaction. This fact provides a link 
to the thermodynamic calculation of Riedel [14], where the correction to the standard 
quasiparticle entropy of a Fermi liquid is presented in the form of the effective bosonic 
contribution. We remind that Riedel [14] summed all the series of ring diagrams rather 
than considered only three first of them as in eq. (4.12). At low temperatures, S^ cm ~ 
T 3 InT [15], i.e. it gives the leading correction to the standard quasiparticle entropy. This 
is the famous correction to the specific heat of liquid 3 He [16,14,15]. Since this correction 
is quite comparable (numerically) to the leading term in the specific heat (~ T), we may 
claim that the liquid 3 He is a liquid with quite strong memory from the point of view of 
kinetics. 

Please also notice that different but an even more simple form of thermodynamic entropy 
including memory corrections follows from eq. (3.46) if one uses thermodynamic relation 
E + PV — fiN = TSV: 

(6.14) 



TS = Tr 



dp 
J2^Y 



An(e) 



6.2 Non-Markovian Entropy in Equilibrium 



We now evaluate the memory correction (s^ em ) eq to the kinetic entropy flow (see eq. 
(5.28)) in thermal equilibrium. Proceeding from the definition of L tj (4.26), as well as 
from equilibrium relations (B.2)-(B.4), and identity (B.5), we present L tj in the form 
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L" + (p) = in h (u)F h (p), L + ~(p) = i [1 + n b (w)] T b (p), 



(6.15) 



where u = p u U u , r b is defined by eq. (6.13), and n b is the bosonic occupation number. 
Now, (Sn lcm ) C q of eq. (5.28) takes the form 



'mcm/cq ■ 



V 

6 



n h (u) [1 + n h (u)\ Flip) 



d A p 
d 4 

^n b (l+n b ) ^-1 
d 4 p 



T 



- 1 



_d_ 

dp v 



[n b (w)r b (p)] (6.16) 



r 



3 9nb + 19T| 
dp,, 3 <9p„ 



n b 



(27T) 



4 ^b 



n b (l + n b ) — 



(6.17) 



The last line is obtained with the help of the partial integration, explicitly taking deriva- 
tives of bosonic occupations and using their property dn b /du; = — n b (l +n b )/T. Now, we 
change the integration variable p — > — p in the last line of eq. (6.17) and use the parity 
properties of the relevant quantities 



F h (-uj) = -F h (uj), n h {-uj) = - [1 + n h (uj)\ , 
and arrive at 



(6.18) 



( S u ) =--u u ( d4p r 3 



(2tt)< 



n b (l + n b ) — 



(6.19) 



Taking the half-sum of the r.h.s. in (6.19) and the last line of eq. (6.17) and using the 
identity (B.7) we finally arrive at 



1 mcmjcq 12 J (2tt) 4 



d 4 p 3 dn h 
(2tt) 4 b <9T' 



(6.20) 



The 0-component of (6.20) precisely coincides with the thermodynamic quantity S°® em 
given by eq. (6.12) in the quasiparticle approximation, provided we consider it in the rest 
frame of the matter, i.e. at U v — {1, 0} . This fact again justifies the name "memory" for 
the thermodynamic quantity S mem . 

Thus, we have demonstrated that our kinetic entropy, including memory contributions, 
coincides with the thermodynamic entropy in thermal equilibrium. 



7 Conclusion and Prospects 



Generalized kinetic equations have been derived from the Kadanoff-Baym equations by 
means of the gradient approximation. Thereby the usual restriction to small mass widths 
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(mass-shell condition) for the particles involved is avoided. The equations are valid, if 
all phase-space distribution functions vary slowly across the space-time region. Since the 
mass-shell relation between energy and momentum is lost for unstable or virtual par- 
ticles, the generalized kinetic equation is formulated in terms of a distribution function 
depending on energy and spatial momentum. Mass-shell equation loses its normal sense in 
context of quasiparticle approximation and rather becomes equivalent to the generalized 
kinetic equation. Alongside a local retarded equation has to be solved which provides the 
dynamical information about the spectral functions of the particles. Besides the usual 
drift and collision terms present in any transport equation, like in Landau's Fermi liquid 
theory [24,4,16], a genuine width and a fluctuation- dependent term appear in the general- 
ized transport equation. The latter term is normally dropped along with the quasiparticle 
approximation which we avoided. This extra term gives rise to backflow contribution in 
the current and produces width and fluctuation dependent contributions to the energy 
momentum tensor pertaining to this transport scheme. It is indeed essential in order to 
preserve the conservation laws in the case of broad damping widths. 

We suggest to follow Baym's ^-derivable principle to consistently construct all terms for 
both, the kinetic and the retarded equation. It has a couple of important advantages, cf. 
ref. [2]. First, it leads to closed, i.e. self-consistent equations which can be closed at any 
order or loop level of the diagrams of <3>, this way defining an effective theory. We showed in 
ref. [8] that the original properties of $ are also valid for genuine nonequilibrium systems 
described within the real-time formalism, namely that the so constructed approximation 
is conserving and at the same time thermodynamically consistent. In fact, for this to 
hold Baym showed that the ^-derivable scheme is not only sufficient but even necessary! 
The ^-derivable energy momentum tensor has explicitly been constructed for the contour 
Dyson equation in ref. [8]. Here we showed that the conserving properties also hold for 
the here derived generalized kinetic equations. As further shown in ref. [8] the scheme 
can be easily extended to include classical boson fields dynamically, such as mean fields 
or condensates. The latter permit to include soft modes in terms of these classical fields, 
much in the spirit of the kinetic picture [34,35,17] of hard loop approximations [36]. 

The structure of the collision term was studied by means of representing the <£> functional 
in terms of " — h" and "H — " Green's functions which represent the Wigner phase-space 
densities [17]. The advantage of the " — h" and "H — " representation is that it leads to a 
natural decomposition of the collision term into multi-particle processes with Feynman 
transition amplitudes which determine the partial rates. Furthermore, it has been dis- 
cussed that sometimes it appears advantageous to account for such memory effects by 
including new "artificial" particles (e.g. bosonization of particle-hole excitations) which 
then lead to local collision terms. 

We also addressed the question whether a closed nonequilibrium system approaches the 
thermodynamic equilibrium during its evolution. Investigating the structure of the colli- 
sion term in the ^-derivable scheme we obtained definite expressions for the local (Marko- 
vian) entropy flow and were able to explicitly demonstrate the if -theorem for some of the 
common choices of $ approximations. The expression for the local entropy flow holds 
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beyond the quasiparticle picture, and thus generalizes the well-known Boltzmann expres- 
sion for the kinetic entropy. To demonstrate memory effects in the generalized kinetics, 
we considered a particular case of a system of fermions interacting via two-body zero- 
range potential. We calculated the memory (non-Markovian) contribution to the kinetic 
entropy, which merges the equilibrium limit with its famous correction to the specific heat 
of liquid 3 He [14-16]: ~ T 3 lnT. Since this correction is quite comparable (numerically) 
to the leading term in the specific heat (~ T), one may claim that the liquid 3 He is a 
liquid with quite strong memory effects from the point of view of kinetics. 

The mass-width effects are important for description of various physical systems. As for 
immediate applications of the developed formalism, we see the description of wide reso- 
nances (such as p-meson, A-resonance, etc.) in nonequilibrium hadron matter produced in 
heavy-ion collisions. Since the widths of these resonances are of the order (or even larger) 
than the excitation energy in the system, a self-consistent treatment of these widths is 
required. Up to now, width effects were considered either within some simplified dynam- 
ics with phenomenological Landau-Migdal residual interaction [37,38] or within a simple 
^-derivable approximation at thermal equilibrium and in the dilute gas limit [31]. In par- 
ticular, it was demonstrated in ref. [17], that the soft-photon production is sensitive to 
dynamical and width effects. Interplay between the width and in-medium population of 
p-mesonic states may also simulate a shift of the p-meson mass in the nuclear medium, 
and thus affects the production of di-leptons in relativistic heavy-ion collisions. It is of 
interest to study these effects within a dynamical approach, such as the scheme presented 
here. 

Further applications concern relativistic plasmas, such as QCD and QED plasmas. The 
plasma of deconfined quarks and gluons was present in the early Universe, it may exist in 
cores of massive neutron stars, and may also be produced in laboratory in ultra-relativistic 
nucleus-nucleus collisions. All these systems need a proper treatment of particle trans- 
port. Perturbative description of soft-quanta propagation suffers of infrared divergences 
and one needs a systematic study of the particle mass-width effects in order to treat 
them, cf. ref. [17]. A thermodynamic ^-derivable approximation for hot relativistic QED 
plasmas — a gas of electrons and positrons in a thermal bath of photons — was recently 
considered by Vanderheyden and Baym [39]. Their treatment may be also applied to 
the high-temperature super-conductors and the fractional quantum Hall effect [40,41]. 
Our approach allows for a natural generalization of such a ^-derivable schemes to the 
dynamical case. 

Another application, as we see it, concerns the description of the neutrino transport in 
supernovas and hot neutron stars during first few minutes of their evolution. At an initial 
stage, neutrinos typically of thermal energy, produced outside (in the mantel) and inside 
the neutron-star core, are trapped within the regions of production. However, coherent 
effects in neutrino production and their rescattering on nucleons [17] reduce the opacity 
of the nuclear-medium and may allow for soft neutrinos to escape the core and contribute 
to the heating off the mantle. The extra energy transport may be sufficient to blow off 
the supernova's mantle in the framework of the shock-reheating mechanism [42]. The 
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description of the neutrinos transport in the semi-transparent region should therefore be 
treated with the due account of mass-widths effects. 
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APPENDICES 

A Matrix Notation 

In calculations that apply the Wigner transformations, it is necessary to decompose the 
full contour into its two branches — the time- ordered and anti-time- ordered branches. One 
then has to distinguish between the physical space-time coordinates x, . . . and the cor- 
responding contour coordinates x c which for a given x take two values x~ = (x~) and 
x+ = ( x ^) (t 1 e {0)1 5 2,3}) on the two branches of the contour (see figure 1). Closed 
real-time contour integrations can then be decomposed as 



where only the time limits are explicitly given. The extra minus sign of the anti-time- 
ordered branch can conveniently be formulated by a { — h} "metric" with the metric 
tensor in { — h} indices 



which provides a proper matrix algebra for multi-point functions on the contour with "co"- 
and "contra" -contour values. Thus, for any two-point function F, the contour values are 
defined as 




(A.l) 




(A.2) 
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F ij (x,y) 
Fi J (x,y) 
Fij(x,y) 



= F(x\y j ), i,j G {-,+}, with 
= a ik F kj (x, y), Fjix, y) := F ik (x, y)a ki 
= a ik a jl F kl (x,y), of = 5., 



ik 



(A.3) 



on the different branches of the contour. Here summation over repeated indices is implied. 
Then contour folding of contour two-point functions, e.g. in Dyson equations, simply 
becomes 

H{x\y k ) = H lk (x,y) = J dz c F(x\z c )G(z c ,y k ) = J dzF^x, z)G jk (z,y) (A.4) 



in the matrix notation. 

For any multi-point function the external point x max , which has the largest physical time, 
can be placed on either branch of the contour without changing the value, since the 
contour-time evolution from a;~ aa , to x+ nax provides unity. Therefore, one-point functions 
have the same value on both sides on the contour. 

Due to the change of operator ordering, genuine multi-point functions are, in general, 
discontinuous, when two contour coordinates become identical. In particular, two-point 
functions like iF(x,y) = (T c A(x)B(y)) become 



\F(x,y) = 



f iF—(x,y) iF- + (x,y)^ 
\iF+-(x,y) iF ++ (x,y)J 



TA{x)B{y) 



T{B{y)A{x) 



\ (A(x)B(y) 



T- l A(x)B(y))J 



, (A.5) 



where T and T _1 are the usual time and anti-time ordering operators. Since there are 
altogether only two possible orderings of the two operators, in fact given by the Wightman 
functions F~ + and F + ~, which are both continuous, not all four components of F are 
independent. Eq. (A.5) implies the following relations between nonequilibrium and usual 
retarded and advanced functions 



F R (x,y) = F—(x,y) - F- + (x,y) = F + ~(x,y) - F ++ (x,y) 

■= <3>(x - y ) (F + -(x, y) - F- + (x, y)) , 
F A (x,y) = F—{x,y) - F + ~(x,y) = F- + (x,y) - F ++ (x,y) 

:= -6(1/0 - x ) (F + -(x, y) - F- + (x, y)) , (A.6) 

where Q(xq — yo) is the step function of the time difference. The rules for the co-contour 
functions F etc. follow from eq. (A.3). 

For such two point functions complex conjugation implies 
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{iF- + (x,y)Y = iF- + (y,x) iF- + (X,p) = real, 

(x, y)) * = (j/, x) =>> iF + ~ {X, p) = real, 
(iF—(x,yj)* = iF ++ (y,x) (iF~(X,p))* = iF++(X,p), 

(F«(x,y))* = F^(y,x) =► (F R (X,p))* = F A (X,p), (A.7) 

where the right parts specify the corresponding properties in the Wigner representation. 
Diagrammatically these rules imply the simultaneous swapping of all + vertices into — 
vertices and vice versa together with reversing the line arrow-sense of all propagator lines 
in the diagram. 

In components the determination of the self-energy from the functional variation of $ (cf. 
eq. (2.9)) reads 



5i$ 5i$ 

the right expression given in the Wigner representation. Note that the variation over 
a {H — }- "contra- variant" Green's function G kl produces a {H — }- "covariant" self-energy 
Sjfc, cf. (A. 3). This occurs due to the same reason discussed above, cf. eq. (A.l), when 
dealing with {H — } matrix notation with integrations over physical times, rather than 
contour times. The extra minus signs occurring for the anti-time ordered branches are 
precisely taken into account by the "covariant" notation (A. 3). 



B Equilibrium Relations 



For completeness of the thermodynamic consideration, we explicitly present here equilib- 
rium relations between quantities on the real-time contour. Basically, they follow from 
the Kubo-Martin-Schwinger condition [43] 



(B.l) 



where e = p u U u — (j, with U u and fj, being a global 4-velocity of the system and a chemical 
potential related to the charge, respectively. All the Green's functions can be expressed 
through the retarded and advanced Green's functions 



G ik {p) = 



[lTn(e)]G R (p)±n(e)G A (p) 



±m(e)A(p) 



\ 



V — i [1 =F n{e)] A(p) - [1 T n(e)\ G A (p) T n(e)G R (p) j 



, (B.2) 
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i, k mean + or — , and the self-energies take a similar form 



£»*(p) 



Here 



/ ^ R (p)±in(e)T(p) =Fira(e)r(p) N 
v i[l T n( £ )]r(p) -Y l A {p)±m{e)T{p) ) 

n(e) = [exp(e/T) ± 



(B-3) 



(B.4) 



are thermal Fermi/Bose-Einstein occupations. They obey some useful relations between 
fermion m and boson n b occupation numbers, like 

n fib (e + u/2) [1 =f nf )b (e - u/2)\ = [n ith (e - u/2) - n ith (e + cu/2)} n h (u), (B.5) 



or derivatives with respect to T 



dn(e) da(e) 



dT 

dn(e) 

dT 



de 



me) 



=F[1 =F n{e)\ ln[l =p - n(e) In 71(e), 



— n b (l=Fn b ). 



(B.6) 
(B.7) 



The link between the Matsubara technique and the real-time formulation used here can be 
provided by extending the real time contour by an imaginary tail going to — i/3, f3 = 1/T, 
this way defining the equilibrium density operator exp(-flH). Thus, the link is pro- 
vided by considering analytic expressions like the "contour trace" of two-point functions 
iF eq (x,y) = (T c A(x)B(y)) 



oo 

CLLO 



J F cq (t,t + 0)dt= J F-+(t,t)dt = ^ J n{e) (F e » - F e ») 

oo 

F Matsubara {ie m ). (B.8) 



C{cq} 



S / 



2ir ie m + //ab - ^ 



Here the contour time t + is placed infinitesimally behind t on the contour in order to 
specify a fixed operator ordering of the two external operators of F, \xab is the chemical 
potential associated to (^A(x) B (y)^, and hab = ~^ba- The step towards the discrete 
Matsubara sum is provided by standard residue technique, cf. ref. [44], fig. 25.4 ff. The 
sum runs over the Matsubara energies 
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(2m + l)nT for fermions 
J V ; (B.9) 

2mnT for bosons. 



1 



Thereby, the Matsubara form of the two-point function F eg has the spectral representation 



7 du F s (u) F ^ z + ^ for lmz>0 

pMatsubara t z \ = / _ cgV 1 " ) _ J f-Q \Q\ 

-°° [ -F C q(-2 + //as) for lmz<0 

in terms of the rea/ iime contour spectral function = — 2ImF C q. 



C Diagram rules 



For relativistic theories with local vertices the diagrammatic rules on the contour are 
identical to the standard Feynman rules except that all time integrations are to be replaced 
by contour integrations. The diagrams contributing to $ are calculated as diagrams with 
one external point, namely the interaction part of the Lagrangian (^C mt (x)^, then contour 
integrated over x and weighted with l/n\, where ri\ counts the number vertices in the 
diagram. These diagrams have to be two-particle irreducible with all lines representing 
full propagators. For details about the corresponding diagrammatic rules on the contour 
see ref. [8]. 

The rules for nonrelativistic two-body interactions are also naturally extended to the 
contour C with 

H hi \ti) = \j<?xij d 4 x 2 ^(x 1 )^(x 2 )V(x 1 - x 2 )^(x 2 )^( Xl ) (C.l) 
c 

with V(xi - x 2 ) = U{xi - x 2 ) Scih - t 2 ) (C.2) 

now defined for contour times t±,t 2 . One has however to observe that for the instantaneous 
two-body interactions the diagrammatic elements are given by <_n_n_r* = — \V(xi — x 2 ), 
i.e. without the factor 1/2 contained in (C.l). With respect to the topological rules, indeed, 
the interaction lines are treated like bosons in the relativistic theory with local vertices. 
To this extend we recall that the expectation value for (^H mt (t)^ — \ x {diagrams} has an 
explicit factor \ in front of the corresponding diagrams. Correspondingly, for the ri\ vertex 
counting for the diagrams of $ each two-body interaction line counts as two vertices, cf. the 
rules given in ref. [4] for the thermodynamic potential in Matsubara formalism. However, 
the two-particle irreducibility of ^-diagrams is defined with respect to cutting dynamical 
propagators only, i.e. leaving the ^-interaction lines untouched. In { — h}-matrix notation 
the above rules imply, cf. [21] sect. X, 
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( ±ns\r* J = *s\ns* = -iV(xi - x 2 ) 



<_n_n_/-t = ~tn_n_r» = 0. (C.3) 
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